clc;
syms x a; 
syms y positive; 
f=exp(-1/(y^2+x^2))*sin(x)^2/x^2*(1+1/y^2)^(x+a^2*y^2); 
L=limit(limit(f,x,1/sqrt(y)),y,inf)